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ABSTRACT 

We study the dynamical evolution of globular clusters using our Henon-type 
Monte Carlo code for stellar dynamics including all relevant physics such as 
two-body relaxation, single and binary stellar evolution, Galactic tidal strip- 
ping, and strong interactions such as physical collisions and binary mediated 
scattering. We compute a large database of several hundred models starting 
from broad ranges of initial conditions guided by observations of young and 
massive star clusters. We show that these initial conditions very naturally 
lead to present day clusters with properties including the central density, core 
radius, half-light radius, and cluster mass, that match well with those of the 
old Galactic globular clusters. In particular, we can naturally reproduce the 
bimodal distribution in observed core radii separating the "core-collapsed" vs 
the "non core-collapsed" clusters. We see that the core-collapsed clusters are 
those that have reached or are about to reach the equilibrium "binary burn- 
ing" phase. The non core-collapsed clusters are still undergoing gravo-thcrmal 
contraction. 

Key words: Galaxy: kinematics and dynamics - Galaxies: star clusters: general 
- globular clusters: general - Methods: numerical 



1 INTRODUCTION 

Studying the evolution of dense star clusters, such as old globular clusters (GCs), is of 
great interest for a variety of branches of astronomy and astrophysics. The high central 



2 S. Chatterjee et al. 

densities and high masses of GCs make them hotbeds for strong dynamical interactions 
facihtating formation of many exotic sources (e.g., X-ray binaries, milhsecond radio pulsars, 
type la supernovae, and blue straggler stars). GCs are important targets for extragalactic 
Astronomy. Detailed observations of their spatial distribution in a galaxy can constrain, for 
example, the potential of the dark matter halo, and give clues to the assembly history of the 
galaxy. The old ages of GCs provide a direct window into the major star-formation episodes 
in the early universe. 

One long-standing question regarding GC s concerns the nature o f their progenitors. The 



observed young rnassive star clusters (e.g., 



1995 



Miller et al 



1997 



Scheepmaker et al. 



Holtzman et al. 



2007 



200^ seem to be potential progenitors 



1992 



Whitmore fc Schweizer 



of the current GCs. The masses, typical sizes, and inferred diss olution timescales 



so-called "super star clusters", as observed, for example, in M51 (I Scheepmaker et al 



or the 



20071 ). 



make them especially good candidates. Interestingly, similar to the Galactic GCs (GGC), 
the sizes of obser yed super star clusters are typ ically a few parsecs independent of the 



cluster mass (e.g., IScheepmaker et al. 



2007 



20091). However, the link between these super 



star clusters and the old evolved GGCs as well as GCs in other galaxies remains speculative. 

The main difficulty in establishing a link between the two populations is that almost 
a Hubble time of evolution separates them. Although the individual qualitative effects of 
each physical proce ss in a GC has been known from decades of numerical studies (e.g.. 



Heggie &: Hutll2003l ). it is impossible to estimate the collective effect of these interdependent 



processes unless a self-consistent simulation is done including them all. For example, two- 
body relaxation leads to a slow energy diffusion from the core to the halo in a GC leading 
to a slow and steady contraction of the core until the gravo-thermal instability occurs and 
the core collapses under gravity. Since, due to equipartition of energy, the low-mass stars are 
scattered to higher velocities relative to their high-mass counterparts, two-body relaxation 
also leads to mass segregation in the cluster. The high-mass stars sink to the core and the 
low-mass stars escape to the halo and preferentially get stripped through the tidal boundary 
of the cluster. Because of this preferential loss of low-mass stars, the stellar MF changes with 
time. The stellar MF, in turn, determines the fraction of low and high-mass st ars in a cluster 



as we ll as the average stellar mass, affecting the mass-segregation timescale ( jGiirkan et al 



20041 ). The contraction of the core via two-body relaxation increases its density. The core 



density determines the interaction rate in the core. These rates affect the binary-single 
(BS) and binary-binary (BB) interaction probabilities at a given time. The BS and BB 



Understanding Core-Collapsed Clusters 3 
interactions in turn generate energy in the core and can support the core stopping further 
contraction, thus, affecting the central density. The BS and BB interactions also change the 
orbital properties of the binaries taking part in these scattering interactions. These changes 
in the binary orbits can alter the evolution pathways that would be taken by a given binary 
which in turn changes, for example, the rate of formation of exotic stellar populations such 
as X-ray binaries and blue straggles stars (BSS). Due to this complexity in the evolution 
of dense star clusters, the only way to learn more about the possible initial properties of 
the observed GCs is to perform numerical simulations including all of the above physical 
processes in tandem with reasonable accuracy for a large enough N. 



The study of GCs has a somewhat long and varied hi story during w 



lich numerical simu- 



lations and observations have complemented each other (IHeggie fc Hut 



20031 ). In particular. 



understanding the physical processes in the cores of GGCs has been of prime interest since 
the evolution of these dense clusters are driven mainly by the energy generation in the core 
and the transport of this energy from the core. The GGCs are observed to show a clear 



bimodal distribution of the core size s that separate the so-called "core-collapsed" c 



from the non core-collapsed clusters (IHarris 



1996 



McLaughlin fc van der Marelll2005l ). The 



usters 



core-collapsed clusters in general show a power-law slope in their density profiles near the 
center. In contras t, the dens ity profiles of the non core-collapsed clusters are well described 



by a King profile (jKinglll966l ) and show a clear fiat part near the center. Theoretical analysis 



based on the current estimated relaxation tim es for the GGCs indicate tha t the majority 



of the GGCs should have had a deep collapse (IDjorgovski 



1993 



Harris 



19961 ). Thus, before 



it was found that all GGCs contain dynamically significant numbers of binaries, theoretical 
studies focused on understanding the process of core-collapse via the balance of outward 
diffusion of energy from the core due to tw o-body relaxation an d post-collapse evolution 
due to dynamical formation of binaries (e.g., iHeggie fc HutI 120031 ). After it was observed in 
the early 1990s that all GGCs contain sufficient number of binaries such that they must 
have been born with substantial primordial binary populations, theoretical studies focused 
on properties of clusters in the "binary-burning" phase in which the outward diffusion of 
energy via two-bod y relaxation is balanced by production of energy vi a djTiamical harden- 



ing of binaries (e.g.. IVesperini fc Chernofi 



1994 



Fregeau &: Rasioll2007l ). It was also realized 



that even a small primord ial binary fraction can support the core from deep collapse for 
more than a Hubble time (iFregeau fc Rasid 120071 ). However, comparison between theoreti- 
cal predictions and observations show that the theoretically predicted core radii during the 
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binary-burning phase for these clusters are at lea st an order of magnitude smaller than 



the observed core radii for the bulk of the GGCs (jVesperini fc Chernofi 



energy sources were propos e d to explaiii these large core sizes (e.g 



19m. Additional 



2008 



Chatterjee et al. 



2008 



Trenti et al. 



Mackey et al 



200 



3 



20071 ). however, these scenarios need rather spe- 



cial conditions and seem unlikely to be satisfied by most of the GGCs. Very basic questions 
remain. Does the bimodal distribution for the GGC core sizes separating the core-collapsed 
and non core-collapsed clusters indicate a physical difference between these clusters? What 
dynamical stage are the core-collapsed and non core-collapsed clusters in today? Can the 
large core sizes of the non core-collapsed (majority of the Galactic population) be explained 
without any special conditions simply as a result of ~ 12Gyr of evolution from realistic 
initial conditions? 

In this study we present the results from a large number of computer simulations (224) 
with initial conditions drawn from a multidimensional grid, spanning all relevant parameter 
ranges as suggested by observations of young star clusters, with the goal of reproducing a 
population of old clusters similar to the GGCs. We use CMC, a Henon-type Monte Carlo 
code including all physical processes such as single and binary stellar evolution, two-body 
relaxation, strong encounters comprising physical collisions and binary-mediated scattering, 
and tidal stripping due to the Galactic tidal field. CMC has been extensively tested and 



the results from CMC s 



(e.g. 



Fregeau fc Rasio 



low excellent agreemerit with those from direct iV- body simulations 



2007 



Chatterjee et al. 



20101 : 



Umbreit et al. 



2OI2I). Our goal is to 



find whether starting from realistic initial conditions (including N, stellar mass function, 
central density, compactness parameter, binary fraction, and cluster siz e) typical of the 



observed young massive star clusters (e.g.. IScheepmaker et al. 



2007 



20091 ) and without any 



special treatment naturally leads to cluster properties similar to the observed old GGCs after 
~ 12 Gyr of evolution. In particular, we focus on understanding the bimodal distribution of 
the observed GGC core radii to identify the dynamical stages of the so called core-collapsed 
and non core-collapsed clusters. 

In Section [2] we briefly explain our code and introduce working definitions for key struc- 
tural parameters of a cluster. In Section [3] we describe the multidimensional grid of initial 
parameters explored in this study. In Section H] we present our key results. In Section [5] we 
show our results identifying the dynamical evolutionary state for the observed core-collapsed 
GGCs. Finally we conclude in Section [61 



2 NUMERICAL METHOD 
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We use our Henon-type Cluster Monte Carlo code (CMC) to numerically model star clusters 
with single and binary stars including all physical processes relevant in globular clusters 
such as two-body relaxation, single and binary stellar evolution, strong encounters including 
physical collisions and binary mediated strong scatt ering encounters. This code has been 



deve 



2003 



oped and rigorously tested over the past decade (I Joshi et al. 



Fregeau &: Rasio 



2007 



Chatter] ee et al 



2010 



Umbreit et al 



2000 



2012 



2001 



Fregeau et al. 



Pattabiraman et al. 



2OI2I ). Using a large grid of ii iitial conditions over the 



young massive clusters (e.g.. 



Scheepmaker et al. 



2007 



range 



of values typical for observed 



2OO9I ) we create over 200 detailed 



star-by-star models and evolve them for 12 Gyr. A handful (about 6) of these models reach 
a very deep collapse phase. For these models the CMC time steps become minuscule and the 
code grinds to a halt. We stop our simulations at that point for these clusters. Note that in 
reality, the deep-collapse phase is halted via formation of the so-called "three-body" binaries 
and the cluster enters into the gravo-thermal oscillation phase. Since in CMC we do not yet 
include the possibility of creating new binaries via three body encounters, we do not address 
this phase at this stage. This however, is not a serious limitation. All clusters that reach 
a very deep collapse stage before t he integration is st opped at 12 Gyr, had zero primordial 



Davis et al 



20081). Here, we only include these clusters 



binaries which is not realistic (e.g., 
in our analysis as limiting cases. Even a small non-zero primordial binary fraction {fb, lowest 
fb = 5% is used in our simulations) can stop the cluster core from collapsing through the 
dynamical hardening of binaries preventing very deep core collapse. At this stage the core size 
remains more or less constant, which is also commonly referred to as "binary-burning" stage. 



A more detailed description of the code, and the various qua. 



stages for a cluster's evolution is presented in IChatterjee et al. 



i tative ly different dynamical 

(bold). 



Since one key goal for this study is to compare the properties of our simulated models 
at about 12 Gyr with the properties of observed GGCs, we have to make sure the same 
parameter definitions are used. In particular there are different definitions for the core radius 
and the cluster size commonly used by theoreticians and observers. The three dimensional 
core radius Vc widely used in A^ -body simulations is a density- weighted measure related 



to the virial radius in the core (iCasertano &: Hut 



19851 ). In contrast, the core radius for 



observed clusters, Tc^obs, is often defined as the distance from the center where the projected 
surface brightness profile drops to half its central value. Alternatively, when observations 



6 S. Chatterjee et al. 

with sufficiently high resolution are available, the core radius can also be based on star-counts 
representing the radius where the projected surface number density is half the corresponding 
central value. If the core radius is calculated using the stellar number density profile, we will 
call it r^Q^jg, while the core radius based on the brightness, or luminosity, density proffie 
is denoted by r^^j^g. In real clusters the density profiles can have large scatter, especially 
close to the centre due to Poisson noise and the presence of bright giants that dominate 
the light but are only few in number. Hence, it is common practice to exclude the light 
from the brightest gi ants when calculating the luminosity density profile for real clusters 



Noyola fc Gebhardtl l2006l ). For the purpose of this study we define 



to reduce noise (e.g., 

r^ohs the core radius obtained from a luminosity density proffie excluding giants with 
> 20 Lq. Since the estimated value of rc^obs is strongly dependent on the estimate of 
the peak luminosity density at the centre of a cluster, values of r^^^^ and r^ohs can differ 
significantly. For real clusters even when star counts are available, only stars above certain 
brightness can be counted due to completeness issues. Hence, we define r^o)^ calculated 
using a number density profile including only stars that are on the main-sequence or on the 



giant 
(e.g., 



jranch with mass es > 0.2 Mq, representing the currently achievable completeness limit 



Leigh et al 



20121). We find that the r^^^ values and the r^^^^ values are not much 



different for our simulated clusters. The three dimensional half-mass radius is defined as 
the radius that includes half of the total mass of the cluster. However, for real clusters this 
quantity is not directly accessible, rather, the two dimensional, projected distance enclosing 
half of the total light of the cluster is calculated and is defined as the half-light radius Vhi- 
We further define rf^^^^s half-light radius for our simulated models by calculating the 

same quantity but excluding bright giants with > 20 L©. We find that the and r'f^^^^g 
values show only minor differences. 



3 INITIAL CONDITIONS 



Our choice of initial conditions is based on observations of super star clusters. All our 
simulated clusters have initial virial radii between r„ = 3 - 4 pc independent of other cluster 
parameters. This range corresponds to initial three-dimensional half mass radii, Vh, ranging 
from 2 - 3pc. These values are in agreement with observations of young, massive clusters 
that indicate that the effective cluster sizes are rather insensitive to the cluster mass, as well 



as metallicity (e.g. lAshman fc Zepj|2001t IScheepmaker et al. 



20071 . 120091 ) and have a median 
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value of ~ 3pc. In addition, observations of old massive LMC clusters, old GCs in NGC 
5128, old clusters in M51, as well as the GGCs indicate that the effective cluster radii show 

119^ 



only a weak correlation with the distance from the galactic center Modge 



1984 



Hesser et al. 



Hwang &: Lee 



20081 ) 



1984 



Mateo 



1987 



van den Bergh et al. 



1991 



Harris 



Scheepmaker et al. 



et al 



2007; 



To restrict the huge parameter space to a certain extent we place all our simulated 
clusters in a circular orbit at a moderate Galactocentric distance of = 8.5 kpc. We avoid 
modeling GCs very close to the Galactic center, where the Galactic field is so strong that 
the tidal stellar loss dominates the cluster's evolution. Due to the assumption of spherical 
symmetry, Monte Carlo codes cannot directly model the tear-drop shaped tidal boundary 
of a star cluster. Instead thes e codes use some pres cription based method (see a detailed 



discussion and calibration in 



Chatter] ee et al.l 120101 ) . If tidal dissolution is the dominant 



driver of the GCs evolution, the approximate method may not be accurate enough. Choosing 
a circular orbit for the simulated clusters is a simplification; however, the results should still 
be valid for eccentric orbits with some effective Galactocentric distance (> 8.5 kpc) (e.g.. 



Baumgardt fc Makind l2003l ). The Galactic tidal field, and consequently the initial tida l 



radius, (r^), for the clusters is calculated following Equation [T] (IBaumgardt fc Makindl2003[ ). 



n 



{ 



1/3 



2/3 
G 1 



(1) 



using a Galactic rotation speed Vq = 220 km s ^. Here the cluster mass is denoted by Md. 
For our set of runs we vary the initial number of stars, A^,;, bet ween 4 — 10 x 10^, encom- 



passing the bulk of the GGCs ( iMcLaughlin &: van der Marel 



20051) . The initial positions and 



velocities are sampled from a King model distribution function with dimensionless potential. 
Wo, in the range 4 — 7.5. We vary the initial binary fraction, /t, between — 0.3. The stellar 



masses o: 



the stars, or primaries in case of a binary, are chosen from the IMF presented in 



Kroupal ( 120011 . their Equations 1 and 2) in the stellar mass range 0.1 — 100 Mq. Secondary 



binary companion masses are sampled from a uniform distribution of mass ratios in the 
range 0.1 M© — mp, where irip is the mass of the primary. The semi-major axes, a, for stellar 
binaries are chosen from a distribution fiat in log within physical limits, namely between 5 x 
the physical contact of the components and the local hard-soft boundary. Although initially 
all binaries in our models are hard at their respective positions, some of these hard binaries 
can become soft during the evolution of the cluster. The cluster contracts under two-body 
relaxation and the velocity dispersion increases making initially hard binaries soft. Moreover, 
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binaries sink to the core due to mass segregation where the velocity dispersion is higher than 
at the initial binary positions. We include these soft binaries in our simulations until they 
are naturally disrupted via strong encounters in the cluster. 



4 BASIC STRUCTURAL PARAMETERS OF OUR MODEL CLUSTERS 



Here we present the evolution of some structural properties of the simulated clusters and 
compare them with the same properties of the observed GGCs. For each of these comparison 
plots, the evolution of a certain cluster property is shown together with the distribution of 
the corresponding observationally derived values for the GGC population. Since all our 
simulated models are at a Galactocentric distance of 8.5 kpc, we also show the observed 
distributions for a subset of GGCs satisfying 7 ^ Rq ^ 10 kpc. Our goal here is to simply 
test whether, starting from observationally motivated initial conditions typical for massive, 
young clusters, the final properties of the cluster ensemble naturally attain ranges of values as 
observed in the GGCs. We do not intend to reproduce the present day distribution for these 
properties since this would require to introduce the distribution of cluster initial conditions 
as another parameter, and, consequently, significantly more simulations which is beyond the 
scope of this study. 

Figure [1] shows the evolution of the total mass of our simulated clusters. The initial 
sharp decrease in the cluster mass (Md) is because of the high mass loss rate via winds and 
compact object formation of the massive stars in the GC. Later on, Md decreases at a slower, 
nearly c onstant, rate, caused by a steady stellar mass loss through the tidal boundary of 



the GC ( iBaumgardt fc Making 



2003 



Davis et al. 



20101 ). The histograms on the right show 



the distribution of observed GGC masses. The observed GGC masses are estimated from 



the absolute visual magnitudes (M^,) given in iHarrid (119961 ) using Equation |2] assuming a 
uniform mass to light ratio M/L = 2 for all clusters. 

^ ^ ]^g-(A/„-4.75)/2. 5+0.30103 ^2) 

Note that this is an approximation. The M/L can vary from cluster to cluster a nd can also 
depend on the evolutionary stage of the cluster in question flAnders et al.ll2009l ). However, 
since here we are only interested in reproducing the cluster mass ranges for most of the GGCs 
and do not try to model a particular cluster, this estimate should be appropriate. Figure [1] 
shows that our model clusters have final masses typical for the bulk of the observed GGCs. 
Since these clusters are modeled with the total N and masses typical of the observed GGCs 
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4 8 12 0.02 



t (Gyr) ^ 

Figure 1. Evolution of the star cluster mass for all simulated models. The solid (black), dotted (red), short-dashed (blue), 
long-dashed (green), and dot-dashed (magenta) lines denote models with primordial = 0, 0.05, 0.1, 0.2, and 0.3, respectively. 
The histograms show the mass distributions of the observed GGCs. The solid histogram is for GGCs with Galactocentric 
distances between 7 and lOkpc. The dashe d histogram is for all observed GGCs. The masses of the observed GGCs are derived 
from their integrated V-band magnitudes llHarrislll99§) using Equation [2] 

the model properties can be directly compared with the overall GGC population without 
any need for scaling. 

Figure [2] shows the evolution of the central density in our simulated models. As a dy- 
namically important GC property, the central density (pc) determines the interaction cross- 
sections for strong scattering inside the core, for example, for BB and BS interactions, and 
stellar collisions. These strong interactions in turn modify the properties of the core, through, 
e.g., binary burning. The pc values sharply decrease during the first ~ 1 Gyr of evolution 
as the high-mass stars lose mass via stellar winds, and compact object formation. Followed 
by the sharp decrease pc increases almost linearly over time during the gravo-thermal core 
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C\J I I I I I I I I u_i I I I 

4 8 12 0.2 0.4 

t (Gyr) P 

Figure 2. Same as Figure[T] but for the central density pc- The pc values for the observed GGCs are also shown in histograms. 
The solid and the dashed histograms are from GGC populations selected as in Figure [T] 



contraction stage. The histograms show the central den sities o: 



convert the bolometric luminosity densities presented in iHarrid (119961 ) and assume M/L = 2. 



the o 



Dserved GGCs. Here we 



The range of central densities for our collection of simulated models is compatible with the 
range for observed GGCs. 

Figure |3] shows the evolution of the density weighted three dimensional core radius r^. 
The Tc for all clusters sharply increases initially for up to about 1 Gyr, because of mass loss 
due to stellar evolution. This mass loss happens mainly at the deepest part of the cluster 
potential since the high-mass stars reside near the cluster center due to mass segregation 
and are affected by mass loss the most. The resulting loss of gravitational binding energy 
expands the cluster core. Once the rate of mass loss goes down, the core contracts via 
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4 8 12 0.4 



t (Gyr) ^ 

Figure 3. Evolution of the three dimensional core radius for all simulated clusters with primordial binary fraction = 0-0.3. 
The line colors and styles have the same meaning as in Figure [l] A few (6) clusters with no primordial binaries reach the deep 
collapse phase within the a Hubble time. We stop integrations for those clusters when this phase is reached. The histograms 
show the observed observed core radii for the GGCs. The solid and dashed histograms show GGCs with Galactocentric distances 
between 7-10 kpc and all GGCs. 

diffusive energy transport from the core to the outside through two-body relaxation. The 
gravo-thermal contraction ceases when this energy flow is balanced by the production of 
energy via dynamical hardening of binaries, the binary-burning phase. In our collection of 
models about 20 clusters show clear binary-burning end stages exhibited by a near constant 
core radius at late stages it > 8 Gyr). The bulk of our models are still contracting at 
tci = 12 Gyr. A few of our models with /{, = do not show the binary-burning stage and 
enter deep-collapse after gravo-thermal contraction. This is because we do not account for 
dynamical creation of binaries. However, as mentioned before, = is a limiting case 
and not representative for observed GCs. Even with fb = 5%, we find that the clusters 
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Figure 4. Same as Figure \3\ but for the ratio of the three dimensional core radius to the three dimensional half mass radius 
(rc/r/j). The rc/r/^ values for the observed GGCs are also shown in histograms. The solid and the dashed histograms are from 
GGC populations selected as in Figure [S] 



that complete the slow contraction phase do not suffer deep collapse, rather reach a steady 
binary-burning energy equilibrium phase. 

Another reliably measured and frequently used parameter reflecting the dynamical state 
of th e evolution of GCs is th e ratio between the core radius Vc and the half-mass radius Vh 

Figure m shows the evolution of Vc/rh for all our clusters. The 



le.g. 



Fregeau fc Rasio 



2001 



range of final r^/ values of the simulated clusters agree well with the observed ones in the 
GGC population, producing values at 12 Gyr close to the peak of the observed distribution. 

We should remind the readers, however, that these Vc and Vh values are not exactly the 
quantities observed directly. A more careful comparison between these structural parameters 
in our models and the observed values in the GGCs requires "observing" the models and 
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80 




Figure 5. Histogram of the ratio of the three dimensional code defined core radius Tc to the observed core radius projected 
core radius ^hs for all simulated clusters. The observed core radii obs are estimated by using the radius from the center 
of the cluster where the projected density is half of the maximum density at the center. The different histograms show the 
distribution of values based on different methods to calculate the density distribution. The thin lines arc for observed core radii 
values calculated using a density distribution including all stars in the cluster. Thick lines are for the same but here the density 
distribution is calculated using a subset of stars satisfying O.OI-L0 < L* < 20 L©. These values ensure that only main-sequence 
stars > 0.2 Mq are included and very bright giants are excluded. The solid (black) and dashed (red) histograms are for observed 
core radii values calculated from the stellar number density, and luminosity density distributions, respectively. 



defining the model quantities such as r^. and as the observers define for a real cluster 
(Section [2]). For example, our model Sl.runl2 has the A^-body defined = 2.3 pc. If 
the same model is observed and the core radius is calculated using the surface number 
density distribution using all main-sequence stars and low-luminosity giant stars, the core 
radius r^^^s* = 2.9 (Table [I]). Figure [5] shows the distribution of values for the ratio Vc^ohs/rc 
using Tc^obs values obtained in 4 different ways. The values for the different definitions of 
Tc^obs and the resulting ratios follow different distributions. Thus, if one wants to estimate 
what the three dimensional A^-body defined value is from the observed value of Tc^obs 
and vice versa, one should be careful about how r^obs for that cluster was calculated. In 
general, r^obs calculated using the number density profile of the cluster show lower levels 
of errors. Moreover, the distributions do not change significantly based on which stars were 
included in the sample to calculate the surface number density profile. The code defined Vc 
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Figure 6. Histogram for the ratio of the three dimensional code defined half-mass radius to the two dimensional observed 
half light radius rf^i ahs- The thin and thick lines have the same meaning as in Figure[5] 

is typically between 1-2 times the r^obs calculated using the number density profile for all 
cluster models in our collection. The luminosity density profiles are subject to a lot more 
noise compared to the number density profiles. As a result the rc/r^obs calculated using the 
luminosity density profiles show a larger spread. In addition, which stars were included in 
calculating the luminosity density profile matters in this calculation significantly. Typically, 
the code-defined Vc values are between 0.25 - 1.5 times and 0.5 - 2.5 times the r^^^^ and 
^c^obs calculated using the luminosity density profile and depending on which stars were 
taken into account for the calculations. 

Similarly, only the half-light radius projected on the sky {vhi) is observed in reality and 
not the half mass radius (r/j). For example, for the simulated cluster Si . runl2 the half light 
radius including all stars is = 4.7 pc. If a luminosity cut-off is used for the same cluster 
the half light radius becomes rf^^^i^s — 4.3 pc. The theoretically calculated three dimensional 
half-mass radius for the same cluster at the same age is Vh = 6.9 pc (Table [T]). Note that this 
difference is not simply due to the projection effect. This difference depends intricately on 
the positions of the bright giant stars that dominate the light but are few in number. Figure 



Understanding Core-Collapsed Clusters 15 

[6] shows the distribution of the ratio Vh/rhi calculated using two separate samples of stars, 
i.e., with and without applying a luminosity cut-off. Typically, the three dimensional half 
mass radius Vh is between 1.5 - 2 times the observed r^y^j^g if the luminosity cut described 



above is used. Note that this range includes the expected geometric factor of y3/2 due 
to projection effect. The spread in values is moderately larger if all stars are included in 
calculating r^i oh&- This is due to the increased statistical fluctuations created by the bright 
and low number of giants with > 20 Lq. 



Hurley 



W e are not the first to po int out this discrepancy in definitions. For instance 
(120071 ) and iTrenti et al.l (120101 ) already found that the A^-body definition of rd can differ 
from an observed Tc/r/j by a factor of a few, and our results confirm that. However, Figures 
[5] and [6] are expected to be useful for observers and theorists modeling GCs since the three 
dimensional Tc or Vh can be be easily estimated using the presented distributions if the 
observed values for these are known. 

We now compare the values of the structural parameters for our simulated models cal- 
culated in a similar way as they are calculated for observed GGCs. Figure [7] shows the pc as 
a function of The observed values for all GGCs and GGCs in the Solar neighborhood 

are also shown in the histograms along the respective axes. The two clusters of points show 
'^w.obs Pc values for clusters with initial = 4 and 3pc, respectively. Our collection of 
models shows very similar half-light radii and central densities compared to those properties 
for the bulk of the GGCs. Note that the final r'f^^i,^ values are strongly dependent on the 
initial for a given Galactocentric distance (r^) as evidenced in the clustering of the points 
in Figure [71 Interestingly, for GCs at tq ~ 8.5 kpc the half-light radii remain close to the 
initial according to our models. Thus, the present observed half-light radii can be used 
as an effective indicator for what the initial was for these clusters. The other parts of the 
observed histograms for r^y^^^g can be easily populated if more simulations are performed 
using smaller initial r„. In addition, using smaller will also lead to smaller final r^jy^j^g 
values for these clusters as they cannot expand further once they fill their tidal radii. 

Figure [8] compares the final values of r^^^g* and pc for all our simulated models with 
the observed core radius distribution for the GGCs. As expected, the smaller the r^ohs^ 
the higher the pc- Our simulated models populate a very similar range in r^obs* values as 
is observed in the bulk of the GGCs, which also applies when other definitions of rc,obs are 
used for the comparison. 

Figure [9] shows the scatter plot for the ?"^obs*/''^/iyobs and pc for our models together 
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Figure 7. The observed half-light radius vs the central density (pc) for all simulated clusters. Histograms on either side 

show the r^jj and pc distributions for the observed GGCs. The ri^i ^bs and pc distributions in our models agree well with 
the bulk of the observed GGC values. The solid and the dashed histograms have the same meaning as in Figure [S] 



with the corresponding histograms for the observed GGCs. The values for this dynamically 
important and dimensionless measure for the compactness of the cluster show excellent 
agreement with the values typically found in all observed GGCs. 

Note that, a long standing puzzle has been the apparent discrepancy between the the- 
oretically predicted Vc/rh values and the values observed for the GGCs. Early numerical 
simulations as well as analytical studies expected that most GGCs are in the binary-burning 
stage. However, the simulated Vc/rh values resulting from binary burning have been found 
to be about an order of magnitude smaller than that for the bulk of the observed popu- 



lation (e.g., IVesperini fc Chernofl 



1994 



Fregeau &: Rasid 120071 ) and it was already known 



that this amount of discrepancy cannot be entirely coming from differences in definitions 
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Figure 8. r^'?"' vs pc- The solid and dashed histograms have the same meaning as in Figure [3] 



(e.g. 



Hurley 



2007 



Trenti et al. 



2010l ). Consequently, additional energy sources in the core 



to expand have been investigated. Several studies proposed different additional energy 



generation mechanisms to exp. 



Chatter] ee et al. 



2008 



Fregeau 



ain the large obseryed Tp /r/, values (e.g 



Trenti et al 



2007 



2008 



Mackey et al. 



20081 ). However, these additional energy 



sources require rather special conditions. For example, high central densities are required to 
create a population of high-mass stars vi a physical collisions tha t can then suffer expedited 
mass loss via compact object formation ( IChatterjee et al.ll2008l ). On the other hand, pres- 



ence of an intermediate mass black hole also cannot be common for the GGCs (iTrenti et al. 



20071 ). Our models, in contrast, are generated without any special assumptions using ob- 
servationally motivated initial conditions, and they naturally create a population of model 
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Figure 9. Ratio of the observed core radius to the observed half light radius '"^q'^s^/'^m obs ^^'^ central density pc- The 
histograms have the same meaning as in Figure [T] 



clusters with properties in excellent agreement with the bulk of the observed GGC properties. 
These results indicate that the progenitors of today's GGCs wer e very similar in proper- 



ties t o the young massive clusters observed, for example, in M 51 (IScheepmaker et al 



2007, 



20091 ). Our results also indicate that the majority of today's GGCs have not yet reached the 
binary-burning, energy-equilibrium stage, and are still contracting under energy transport 
via two-body relaxation. 



5 WHAT IS A "CORE-COLLAPSED" CLUSTER? 

A lot of early theoretical work was devoted to understand the gravo-thermal collapse and 
subsequent evolution of the core and the cluster as a whole due to hardening of primordial 
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Figure 10. Examples for a binary-burning cluster (left) and a core-contracting cluster (right) from our models chosen randomly. 
The top panels show the evolution of rc/r^ for each clusters. The bottom panels show the surface luminosity density profiles 
calculated excluding the bright (L^ > 20 Lq) giants to reduce noise. The errorbars are estimated la Poisson errors for each bin. 
The binary-burning model shows a clear power-law slope until the data is too noisy. In contrast, the core-contracting model 
shows a clear King density profile. 



binaries. However, analytical results as well as numerical simulations showed that the Tc 
values in the binary-burning phase are about an order of magnitude too small compared to 
the bulk of the GGC cores. 

Our simulation results presented in Section H] show that this is simply because the bulk 
of the observed GGCs are not in binary-burning equilibrium stage. Rather they are still con- 
tracting. Now we focus on understanding the clearly bimodal distribution of the core radii 
of the GGCs. Depending on the shape of the observed cluster density profiles, all GGCs 
are divided i nto two categories, namely core-col. 



Harris 



1996 



apsed and non core-collapsed clusters (e.g., 
McLaughlin fc van der Marelll2005l ). The so called core-collapsed clusters ex- 



hibit a power-law increase in the density profile until the limit of resolution of the observation. 



The non core-collapsed clusters show a very clear flattenin g of t 



le den sity profile and the 



19661 ). Is there a distinct 



profile for these clusters are well fitted with a King profile fiKingI 
difference dynamically between the two categories of clusters? 

Figure [10] shows examples of two representative clusters chosen randomly from our large 
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Figure 11. Distribution of r'^'^^^/rf'Y* . . The thick and thin lines denote values from our simulated models and observed GGCs, 

^ CjObs / hi, ohs ' 

respectively. The solid and dashed histograms show distributions for model clusters in binary-burning energy equilibrium and in 
gravo-thermal contraction stages, respectively. The solid and dashed histograms for the observed GGCs denote core-collapsed 
clusters and non core-collapsed clusters, respectively. The ratio between the core and half-light radius for binary-burning clusters 
have very similar values compared to the values for the core-collapsed population among the GGCs. On the other hand the 
bulk of the GGCs have values for this ratio similar to the model clusters that are still in the gravo-thermal contraction stage 
at integration stopping time (fs 12Gyr). 



collection of models, one is in a clearly binary-burning stage and the other is still contracting. 
The evolution of rc/r^ is shown for the two clusters as well as their respective surface 
density profiles. The binary-burning cluster can be clearly identified by the near constant 
value of Vc/rfi starting at about 8 Gyr. The core contracting cluster shows a constant rate of 
contraction until 12 Gyr. The surface density profiles for the two clusters are very different. 
The surface density profile for the binary-burning cluster looks very much like a so called 
core-collapsed cluster and exhibits a clearly power-law slope to very small radius below which 
the profile is noisy. In contrast, the surface density profile for the core contracting cluster 
shows a clear King profile with a distinct fiat part near the center below a few parsecs. 

We divide the full sample of our simulations in two subgroups: 1) binary-burning: models 
showing a clear binary-burning stage before integration is stopped, and 2) core-contracting: 
models still contracting due to two-body relaxation until integration stopping time. Figure 
[TT] shows the distributions for r^o^s / ^u^hs ^^e two subgroups of our theoretical models. 
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The distributions peak at different values, with the binary-burning clusters having much 
lower core radii. The distributions for r^obs for the observed core-collapsed and non core- 
collapsed GGCs are also shown. Note that the values for the binary-burning clusters from 
our models and the core-collapsed clusters from the observed GGCs show a very similar 
range of values. Similarly, r^obs / ''^hi^bs the contracting clusters in our models show very 
similar values compared to the those of the observed non core-collapsed GGCs. Also note 
that most of the core-contracting models with small r^^^ < 0.1 pc values are about to 
start binary-burning, but was classified by contracting since the constant Vc/rh stage of 
evolution was not clearly seen. The observed GGCs show larger ranges for the core radii 
values compared for both categories of clusters compared to our models. This is simply 
because we are forced to limit the ranges of the grid of initial conditions to constrain the 
number of required cluster calculations to a tractable amount. We remind the readers that 
the heights of the histograms between the model clusters and the observed GGCs are not 
compared here, since that depends directly on the distributions of initial cluster parameters, 
determination of which is beyond the scope of this study, and are chosen arbitrarily. Only 
ranges in values are compared. 



6 CONCLUSIONS 



We have presented a large (~ 200) collection of cluster models created with the Northwest- 
ern group's Henon-type Monte Carlo code in star- by-star detail. We start our models with 
initial parameters including the stellar mass spectrum, cluster size, concentration, and pri- 
mordial binar y fraction fh over large ranges (Section [3]) guided by observed young massive 



clusters (e.g., IScheepmaker et al. 



2007 



20091 ). We find that these initial clusters very nat- 



urally produce a population of model clusters with structural properties including cluster 
mass, Pc Tc, and Vh in excellent agreement with the bulk of the GGC properties after about 
12 Gyr of evolution without any special considerations or fine tuning (e.g., very high density 



to aid collisional expedited stellar mass loss via compact o bject formation; IChatterjee et al. 



2008 



or intermediate mass black holes; 



Trenti et al. 



20071 ). We pay attention to the various 



different commonly used definitions of the structural parameters Vc and r/j and calculate 
these quantities from our models as an observer would for real clusters. These parameters 
are then compared and found to agree well with the ranges from observed GGCs. Using our 
large collection of models we also show the distribution of the ratio of the three dimensional 
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code-defined Tc and to the corresponding "observed" values. We expect that these distri- 
butions of ratios for the and r/j values will be valuable for observers and theorists alike to 
convert the values of these parameters from one set of definitions to another. 

From the evolution of the code-defined three-dimensional structural parameters of all our 
models, we find that all qualitatively different evolutionary stages are observed, in particular, 
the initial expansion due to stellar evolution driven mass loss, core contraction driven by 
two-body relaxation, and the binary-burning equilibrium stage (for clusters with > 0) 
driven by a balance between energy production via dynamical hardening of binaries in the 
core and outward diffusion of energy from the core due to two-body relaxation. Our results 
indicate that the progenitors of today's GGCs were very sin ailar in properties to the presen t 



day young massive clusters (observed, for example, in M 51 IScheepmaker et al. 



2007 



20091). 



Of course, the metallicities of these progenitors must have been much lower compared to 
today's massive young clusters. 

After establishing that our collection of cluster models are representative of the observed 
GGCs we investigate the apparently bimodal distribution of the observed core radii of the 
GGCs created by the core collapsed clusters and non core-collapsed clusters. In particular, 
we answer the question if the core collapsed clusters are dynamically different from the non 
core collapsed clusters. We find that the surface brightness profile for the binary-burning 
cluster shows a prominent power-law slope near the center. The core collapsed clusters 
are observationally defin ed as the clusters that show this d istinct feature in their surface 



brightness profiles (e.g., iMcLaughlin fc van der Marel 



20051 ) . In contrast, a model cluster 



that is still contracting at 12 Gyr shows a s urface brigh tness profile that has a clear flat 
central part and is fitted well by a King profile (jKinglll966l ). We further divide our models into 
two subsets, one containing clusters in the binary-burning stage, and the other containing 
clusters that are still contracting at 12 Gyr. We compare the ratio between the core and half 
light radii of our binary-burning and contracting clusters with those for the core-collapsed 
and non core-collapsed clusters in the observed GGCs, respectively. We find that the binary- 
burning ?"^obsV''^ft,y,obs values in our models are in agreement with those of the core-collapsed 
GGCs. Similarly, the contracting ?"^obsV'"/i"obs values in our models are in agreement with 
those of the non core-collapsed GGCs (Figure fTTl) . In addition, the surface density profiles for 
the binary-burning and contracting models show features (power-law slope and flat central 
part, respectively) similar to those seen in the surface density profiles of the core-collapsed 
and non core-collapsed GGCs (Figure [T0|) . Thus our results clearly indicate that the so called 
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core-collapsed GGCs are in fact at the binary-burning stage whereas, the non core-collapsed 
GGCs are still contracting under two-body relaxation. This also indicates that the majority 
of the GGCs (since most GGCs are non core-collapsed) are not in energy equ ilibrium as 
was expected by some earlier theoretical models (e.g., iFregeau fc Rasid 120071 ). One key 
implication for this finding is that analytical estimates of interaction rates in a GGC must 
take into account the fact that the present day observed structural parameters including 
Pc has not been constant and is still evolving. Hence, to calculate a correct estimate one 
must integrate the time dependent cross-sections (based on the changi ng values of th ese 
parameters) over an appropriate length of time as has been done in, e.g., iFregeaul (120081 ). 
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from UF Department of Astronomy and College of Liberal Arts and Sciences. 



Table 1: List of simulations: Wq is the centra. 



concen- 



tration parameter for a King profile (jKing]|l966l ). cluster 
mass M is in 10^ M©, number of bound cluster objects 
is given in 10^, central stellar mass density pc is in 
lO^Mo/pc^ and r/j are in pc, time of final snapshot 
t is in Gyr and c is the concentration parameter defined 
as log(rt/rc). All final values are extracted from the final 
snapshot of the simulated clusters. The core radii and 
half-light radii calculated as they are defined in observed 
clusters, r^^hs ^hi^ohs ^Iso listed. 



Name 










Initial 


















Final 












Wo 


M 




rc 


rh 


Pc 


h 


fb,C 


c 


t 


M 


N 


rc 


N,cut 
c,obs 




„cut 
hl,ohs 


Pc 


h 




c 


runl 
run2 
runS 
run4 
run5 


4 
4 
4 
4 
4 


2.5 
3.8 
5.1 
6.4 
2.6 


4 
6 
8 

10 
4 


1.6 
1.6 
1.6 
1.6 
1.6 


3.3 
3.3 
3.3 
3.3 
3.3 


12.2 
16.2 
23.3 
29.0 
12.3 


0.00 
0.00 
0.00 
0.00 
0.05 


0.00 
0.00 
0.00 
0.00 
0.05 


1.3 
1.6 
1.6 
1.7 
1.3 


12 
12 
12 
12 
12 


1.4 
2.1 
2.8 
3.5 
1.4 


3 
5 
7 
9 
3 


2.0 
2.3 
2.4 
2.5 
2.1 


1.7 
1.9 
1.9 
2.9 
2.3 


7.1 
7.0 
6.9 
6.8 
7.2 


4.2 
4.3 
4.4 
4.4 
4.2 


1.0 
1.1 
1.3 
1.5 
0.9 


0.00 
0.00 
0.00 
0.00 
0.05 


0.00 
0.00 
0.00 
0.00 
0.07 


1.5 
1.5 
1.5 
1.5 
1.4 
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Table 1 - Continued 



Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



run6 


4 


3.9 


6 


1.6 


3.3 


17.3 


0.05 


0.05 


1 


4 


12 


2.1 


5 


2.3 


2.1 


7.1 


4.4 


1.1 


0.05 


0.07 


1.5 


run? 


4 


5.3 


8 


1.6 


3.3 


24.7 


0.05 


0.05 


1 


5 


12 


2.9 


7 


2.4 


2.2 


7.0 


4.4 


1.2 


0.05 


0.07 


1.5 


runS 


4 


6.6 


10 


1.6 


3.3 


30.2 


0.05 


0.05 


1 


8 


12 


3.6 


9 


2.5 


3.3 


6.9 


4.4 


1.3 


0.05 


0.06 


1.5 


run9 


4 


2.7 


4 


1.6 


3.3 


12.4 


0.10 


0.10 


1 


5 


12 


1.4 


3 


2.2 


2.1 


7.3 


4.4 


0.8 


0.09 


0.14 


1.4 


run 10 


4 


4.0 


6 


1.6 


3.3 


17.7 


0.10 


0.10 


1 


6 


12 


2.2 


5 


2.3 


1.9 


7.2 


4.4 


1.1 


0.09 


0.14 


1.5 


run 11 


4 


5.4 


8 


1.6 


3.3 


25.3 


0.10 


0.10 


1 


5 


12 


2.9 


7 


2.5 


2.0 


7.1 


4.4 


1.2 


0.09 


0.13 


1.5 


run 12 


4 


6.8 


10 


1.6 


3.3 


30.9 


0.10 


0.10 


1 


5 


12 


3.7 


9 


2.6 


2.0 


7.0 


4.4 


1.3 


0.09 


0.13 


1.5 


runl3 


4 


2.8 


4 


1.6 


3.3 


14.3 


0.20 


0.20 


1 





12 


1.5 


3 


2.3 


2.5 


7.5 


4.5 


0.8 


0.18 


0.26 


1.4 


run 14 


4 


4.3 


6 


1.6 


3.3 


18.3 


0.20 


0.20 


1 


5 


12 


2.3 


5 


2.4 


2.1 


7.3 


4.3 


1.0 


0.18 


0.26 


1.5 


run 15 


4 


5.7 


8 


1.6 


3.3 


27.6 


0.20 


0.20 


1 


7 


12 


3.0 


7 


2.6 


2.2 


7.2 


4.5 


1.1 


0.18 


0.25 


1.5 


run 16 


4 


7.1 


10 


1.6 


3.3 


33.4 


0.20 


0.20 


1 


5 


12 


3.8 


9 


2.7 


2.1 


7.2 


4.5 


1.2 


0.18 


0.24 


1.5 


run 17 


5 


2.5 


4 


1.5 


3.3 


14.2 


0.00 


0.00 


1 





12 


1.4 


3 


1.9 


1.2 


7.2 


4.3 


1.4 


0.00 


0.00 


1.5 


run 18 


5 


3.8 


6 


1.5 


3.3 


18.8 


0.00 


0.00 


1 


8 


12 


2.1 


5 


2.1 


1.7 


7.0 


4.4 


1.4 


0.00 


0.00 


1.5 


run 19 


5 


5.1 


8 


1.5 


3.3 


27.2 


0.00 


0.00 


1 


4 


12 


2.8 


7 


2.3 


1.5 


6.9 


4.4 


1.5 


0.00 


0.00 


1.5 


run20 


5 


6.4 


10 


1.5 


3.3 


33.3 


0.00 


0.00 


1 


6 


12 


3.5 


9 


2.3 


2.4 


6.9 


4.4 


1.7 


0.00 


0.00 


1.5 
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Table 1 - Continued 



Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



run21 


5 


2.6 


4 


1.5 


3.3 


14.3 


0.05 


0.05 


1 





12 


1.4 


3 


2.0 


1.4 


7.3 


4.4 


1.2 


0.05 


0.08 


1.5 


run22 


5 


3.9 


6 


1.5 


3.3 


20.1 


0.05 


0.05 


1 


6 


12 


2.1 


5 


2.2 


1.4 


7.1 


4.4 


1.2 


0.05 


0.07 


1.5 


run23 


5 


5.3 


8 


1.5 


3.3 


28.9 


0.05 


0.05 


1 


3 


12 


2.9 


7 


2.3 


2.0 


7.1 


4.4 


1.4 


0.05 


0.07 


1.5 


run24 


5 


6.6 


10 


1.5 


3.3 


34.7 


0.05 


0.05 


1 


5 


12 


3.6 


9 


2.4 


2.6 


7.0 


4.4 


1.7 


0.05 


0.07 


1.5 


run25 


5 


2.7 


4 


1.5 


3.3 


14.4 


0.10 


0.10 


1 





12 


1.4 


3 


2.0 


1.3 


7.4 


4.4 


1.1 


0.09 


0.14 


1.5 


run26 


5 


4.0 


6 


1.5 


3.3 


20.5 


0.10 


0.10 


1 


2 


12 


2.2 


5 


2.2 


1.4 


7.3 


4.4 


1.2 


0.09 


0.14 


1.5 


run27 


5 


5.4 


8 


1.5 


3.3 


29.6 


0.10 


0.10 


1 


2 


12 


2.9 


7 


2.4 


2.0 


7.2 


4.3 


1.4 


0.09 


0.13 


1.5 


run28 


5 


6.8 


10 


1.5 


3.3 


35.4 


0.10 


0.10 


1 


3 


12 


3.7 


9 


2.4 


2.4 


7.1 


4.5 


1.6 


0.09 


0.13 


1.5 


run29 


5 


2.8 


4 


1.5 


3.3 


16.8 


0.20 


0.20 


1 


2 


12 


1.5 


3 


2.2 


2.1 


7.6 


4.4 


0.9 


0.18 


0.26 


1.4 


runSO 


5 


4.3 


6 


1.5 


3.3 


21.2 


0.20 


0.20 


1 


2 


12 


2.3 


5 


2.3 


1.9 


7.4 


4.5 


1.1 


0.18 


0.26 


1.5 


runSl 


5 


5.7 


8 


1.5 


3.3 


32.3 


0.20 


0.20 


1 


4 


12 


3.0 


7 


2.5 


2.2 


7.3 


4.4 


1.2 


0.18 


0.24 


1.5 


run32 


5 


7.1 


10 


1.5 


3.3 


38.4 


0.20 


0.20 


1 


6 


12 


3.8 


9 


2.5 


2.0 


7.2 


4.4 


1.6 


0.18 


0.25 


1.5 


run33 


6 


2.5 


4 


1.4 


3.2 


17.2 


0.00 


0.00 


1 


3 


12 


1.4 


3 


1.7 


0.9 


7.3 


4.3 


1.8 


0.00 


0.00 


1.5 


run34 


6 


3.8 


6 


1.4 


3.3 


22.4 


0.00 


0.00 


1 


8 


12 


2.1 


5 


1.9 


1.2 


7.1 


4.3 


1.9 


0.00 


0.00 


1.6 


run35 


6 


5.1 


8 


1.4 


3.3 


32.5 


0.00 


0.00 


1 


4 


12 


2.8 


7 


2.1 


1.6 


7.0 


4.3 


2.0 


0.00 


0.00 


1.5 
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Table 1 - Continued 



Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



run36 


6 


6.4 


10 


1.4 


3.3 


40.0 


0.00 


0.00 


1 


5 


12 


3.5 


9 


2.2 


2.1 


7.0 


4.4 


2.2 


0.00 


0.00 


1.6 


run37 


6 


2.6 


4 


1.4 


3.2 


17.3 


0.05 


0.05 


1 


2 


12 


1.4 


3 


1.8 


1.5 


7.4 


4.4 


1.5 


0.05 


0.08 


1.5 


run38 


6 


3.9 


6 


1.4 


3.3 


24.0 


0.05 


0.05 


1 


2 


12 


2.1 


5 


2.0 


1.7 


7.3 


4.4 


1.6 


0.05 


0.07 


1.5 


run39 


6 


5.3 


8 


1.4 


3.3 


34.6 


0.05 


0.05 


1 


2 


12 


2.8 


7 


2.2 


1.7 


7.2 


4.4 


1.7 


0.05 


0.07 


1.5 


run40 


6 


6.6 


10 


1.4 


3.3 


41.6 


0.05 


0.05 


1 


5 


12 


3.6 


9 


2.2 


1.7 


7.1 


4.4 


2.1 


0.05 


0.07 


1.5 


run41 


6 


2.7 


4 


1.4 


3.2 


17.4 


0.10 


0.10 


1 


2 


12 


1.4 


3 


1.9 


1.5 


7.5 


4.4 


1.5 


0.09 


0.15 


1.5 


run42 


6 


4.0 


6 


1.4 


3.3 


24.5 


0.10 


0.10 


1 


2 


12 


2.2 


5 


2.1 


1.3 


7.4 


4.4 


1.5 


0.09 


0.14 


1.5 


run43 


6 


5.4 


8 


1.4 


3.3 


35.4 


0.10 


0.10 


1 


2 


12 


2.9 


7 


2.2 


2.0 


7.3 


4.5 


1.7 


0.09 


0.13 


1.5 


run44 


6 


6.8 


10 


1.4 


3.3 


42.5 


0.10 


0.10 


1 


8 


12 


3.7 


9 


2.3 


1.9 


7.2 


4.5 


1.9 


0.09 


0.13 


1.6 


run45 


6 


2.8 


4 


1.4 


3.2 


20.3 


0.20 


0.20 


1 


2 


12 


1.5 


3 


2.1 


1.1 


7.8 


4.5 


1.1 


0.18 


0.27 


1.4 


run46 


6 


4.3 


6 


1.4 


3.3 


25.2 


0.20 


0.20 


1 


1 


12 


2.3 


5 


2.2 


1.7 


7.6 


4.5 


1.3 


0.18 


0.26 


1.5 


run47 


6 


5.7 


8 


1.4 


3.3 


38.7 


0.20 


0.20 


1 


5 


12 


3.0 


7 


2.4 


1.8 


7.4 


4.4 


1.4 


0.18 


0.24 


1.5 


run48 


6 


7.1 


10 


1.4 


3.3 


46.0 


0.20 


0.20 


1 


5 


12 


3.8 


9 


2.4 


2.4 


7.3 


4.5 


1.8 


0.18 


0.25 


1.5 


run49 


7 


2.5 


4 


1.3 


3.2 


21.6 


0.00 


0.00 


1 


2 


12 


1.4 


3 


1.5 


1.0 


7.4 


4.4 


3.0 


0.00 


0.00 


1.6 


run50 


7 


3.8 


6 


1.3 


3.2 


27.6 


0.00 


0.00 


1 


5 


12 


2.1 


5 


1.7 


1.6 


7.3 


4.4 


2.6 


0.00 


0.00 


1.6 
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Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



run51 


7 


5.1 


8 


1.3 


3.3 


40.2 


0.00 


0.00 


1 


4 


12 


2.8 


7 


1 


8 


1.5 


7.2 


4.4 


3.0 


0.00 


0.00 


1.6 


run52 


7 


6.4 


10 


1.3 


3.2 


49.6 


0.00 


0.00 


1 


6 


12 


3.5 


9 


1 


9 


1.2 


7.1 


4.4 


3.1 


0.00 


0.00 


1.6 


run53 


7 


2.6 


4 


1.3 


3.2 


21.7 


0.05 


0.05 


1 


4 


12 


1.4 


3 


1 


5 


1.2 


7.6 


4.5 


2.5 


0.05 


0.08 


1.6 


run54 


7 


3.9 


6 


1.3 


3.2 


29.7 


0.05 


0.05 


1 


5 


12 


2.1 


5 


1 


8 


1.5 


7.4 


4.5 


2.4 


0.05 


0.08 


1.6 


run55 


7 


5.3 


8 


1.3 


3.3 


43.0 


0.05 


0.05 


1 


2 


12 


2.8 


7 


1 


9 


1.4 


7.3 


4.4 


2.6 


0.05 


0.07 


1.6 


run56 


7 


6.6 


10 


1.3 


3.2 


51.5 


0.05 


0.05 


1 


6 


12 


3.6 


9 


2 





1.7 


7.2 


4.5 


2.9 


0.05 


0.07 


1.6 


run57 


7 


2.7 


4 


1.3 


3.2 


21.7 


0.10 


0.10 


1 


2 


12 


1.4 


3 


1 


7 


1.3 


7.7 


4.5 


1.8 


0.09 


0.16 


1.5 


run58 


7 


4.0 


6 


1.3 


3.2 


30.3 


0.10 


0.10 


1 


3 


12 


2.2 


5 


1 


9 


1.5 


7.5 


4.5 


2.0 


0.09 


0.14 


1.5 


run59 


7 


5.4 


8 


1.3 


3.3 


44.0 


0.10 


0.10 


1 


4 


12 


2.9 


7 


2 





1.4 


7.4 


4.5 


2.5 


0.09 


0.14 


1.6 


run60 


7 


6.8 


10 


1.3 


3.2 


52.5 


0.10 


0.10 


1 


4 


12 


3.6 


9 


2 


1 


1.4 


7.3 


4.5 


2.5 


0.09 


0.13 


1.6 


run61 


7 


2.8 


4 


1.3 


3.2 


25.5 


0.20 


0.20 


1 


2 


12 


1.5 


3 


1 


9 


1.4 


7.9 


4.6 


1.4 


0.18 


0.27 


1.5 


run62 


7 


4.3 


6 


1.3 


3.2 


31.2 


0.20 


0.20 


1 


4 


12 


2.3 


5 


2 


1 


1.5 


7.7 


4.5 


1.6 


0.18 


0.26 


1.5 


run63 


7 


5.7 


8 


1.3 


3.3 


48.2 


0.20 


0.20 


1 


3 


12 


3.0 


7 


2 


1 


1.9 


7.6 


4.5 


2.1 


0.18 


0.26 


1.5 


run64 


7 


7.1 


10 


1.3 


3.2 


56.9 


0.20 


0.20 


1 


7 


12 


3.8 


9 


2 


2 


1.7 


7.5 


4.5 


2.1 


0.18 


0.24 


1.5 


run65 


8 


2.5 


4 


1.2 


3.2 


28.5 


0.00 


0.00 


1 


6 


12 


1.4 


3 


1 


1 


0.8 


7.7 


4.5 


7.1 


0.00 


0.00 


1.7 
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Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



run66 


8 


3.8 


6 


1.2 


3.2 


35.7 


0.00 


0.00 


1 


3 


12 


2.1 


5 


1.4 


0.8 


7.5 


4.5 


4.7 


0.00 


0.00 


1.7 


run67 


8 


5.1 


8 


1.2 


3.2 


52.8 


0.00 


0.00 


1 


8 


12 


2.8 


7 


1.6 


1.2 


7.4 


4.5 


4.7 


0.00 


0.00 


1.7 


run68 


8 


6.4 


10 


1.2 


3.2 


64.3 


0.00 


0.00 


1 


8 


12 


3.5 


9 


1.6 


1.2 


7.3 


4.5 


5.1 


0.00 


0.00 


1.7 


run69 


8 


2.6 


4 


1.2 


3.2 


28.6 


0.05 


0.05 


1 


4 


12 


1.4 


3 


1.3 


1.0 


7.8 


4.5 


4.0 


0.05 


0.08 


1.6 


run70 


8 


3.9 


6 


1.2 


3.2 


38.5 


0.05 


0.05 


1 


3 


12 


2.1 


5 


1.5 


1.1 


7.7 


4.6 


4.1 


0.05 


0.08 


1.6 


run71 


8 


5.3 


8 


1.2 


3.2 


56.6 


0.05 


0.05 


1 


5 


12 


2.8 


7 


1.6 


1.3 


7.5 


4.6 


4.1 


0.05 


0.08 


1.6 


run72 


8 


6.6 


10 


1.2 


3.2 


66.5 


0.05 


0.05 


1 


7 


12 


3.6 


9 


1.8 


1.4 


7.4 


4.5 


4.2 


0.05 


0.07 


1.6 


run73 


8 


2.7 


4 


1.2 


3.2 


28.6 


0.10 


0.10 


1 


2 


12 


1.4 


3 


1.4 


1.1 


7.9 


4.6 


3.4 


0.09 


0.16 


1.6 


run74 


8 


4.0 


6 


1.2 


3.2 


39.3 


0.10 


0.10 


1 


3 


12 


2.2 


5 


1.6 


1.0 


7.8 


4.6 


3.2 


0.09 


0.15 


1.6 


run75 


8 


5.4 


8 


1.2 


3.2 


57.9 


0.10 


0.10 


1 


4 


12 


2.9 


7 


1.8 


1.5 


7.6 


4.6 


3.3 


0.09 


0.14 


1.6 


run76 


8 


6.8 


10 


1.2 


3.2 


67.9 


0.10 


0.10 


1 


8 


12 


3.6 


9 


1.8 


1.4 


7.5 


4.6 


3.9 


0.09 


0.14 


1.6 


run77 


8 


2.8 


4 


1.2 


3.2 


33.7 


0.20 


0.20 


1 


2 


12 


1.5 


3 


1.7 


1.2 


8.2 


4.6 


2.1 


0.18 


0.28 


1.5 


run78 


8 


4.3 


6 


1.2 


3.2 


40.4 


0.20 


0.20 


1 


2 


12 


2.3 


5 


1.8 


1.5 


8.0 


4.7 


2.7 


0.18 


0.27 


1.6 


run79 


8 


5.7 


8 


1.2 


3.2 


63.4 


0.20 


0.20 


1 


7 


12 


3.0 


7 


1.9 


1.5 


7.8 


4.6 


2.9 


0.18 


0.26 


1.6 


runSO 


8 


7.1 


10 


1.2 


3.2 


73.6 


0.20 


0.20 


1 


7 


12 


3.8 


9 


2.0 


1.4 


7.7 


4.7 


3.2 


0.18 


0.25 


1.6 
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Name Initial 




Final 






Wo M N Th Pc 


fb fb, c c 


t M N Tc <r 


Pc 


fb fb, c c 



runSl 


6.5 


2.5 


4 


1.1 


3.2 


40.1 


0.00 


0.00 


1.4 


12 


1.4 


3 


0.5 


0.7 


8.0 


4.7 


64.9 


0.00 


0.00 


2.1 


run82 


6.5 


3.8 


6 


1.1 


3.2 


49.3 


0.00 


0.00 


1.5 


12 


2.1 


5 


1.0 


0.6 


7.8 


4.6 


14.7 


0.00 


0.00 


1.8 


run83 


6.5 


5.1 


8 


1.1 


3.2 


73.0 


0.00 


0.00 


1.6 


12 


2.8 


7 


1.2 


0.8 


7.6 


4.7 


11.7 


0.00 


0.00 


1.8 


run84 


6.5 


6.4 


10 


1.1 


3.2 


88.8 


0.00 


0.00 


2.0 


12 


3.5 


9 


1.3 


0.8 


7.5 


4.6 


10.0 


0.00 


0.00 


1.8 


run85 


6.5 


2.6 


4 


1.1 


3.2 


40.3 


0.05 


0.05 


1.6 


12 


1.4 


3 


0.9 


0.5 


8.1 


4.7 


14.0 


0.05 


0.10 


1.8 


run86 


6.5 


3.9 


6 


1.1 


3.2 


53.4 


0.05 


0.05 


1.4 


12 


2.1 


5 


1.2 


0.8 


7.9 


4.7 


8.4 


0.05 


0.08 


1.7 


run87 


6.5 


5.3 


8 


1.1 


3.2 


78.5 


0.05 


0.05 


1.5 


12 


2.8 


7 


1.3 


0.8 


7.8 


4.8 


9.0 


0.05 


0.08 


1.7 


run88 


6.5 


6.6 


10 


1.1 


3.2 


91.7 


0.05 


0.05 


1.6 


12 


3.6 


9 


1.4 


1.0 


7.7 


4.7 


8.7 


0.05 


0.08 


1.7 


run89 


6.5 


2.7 


4 


1.1 


3.2 


40.2 


0.10 


0.10 


1.4 


12 


1.4 


3 


1.1 


0.8 


8.2 


4.7 


6.1 


0.09 


0.17 


1.7 


run90 


6.5 


4.0 


6 


1.1 


3.2 


54.4 


0.10 


0.10 


1.4 


12 


2.2 


5 


1.3 


0.9 


8.0 


4.7 


5.8 


0.09 


0.16 


1.7 


run91 


6.5 


5.4 


8 


1.1 


3.2 


80.2 


0.10 


0.10 


1.6 


12 


2.9 


7 


1.4 


0.9 


7.9 


4.8 


7.3 


0.09 


0.15 


1.7 


run92 


6.5 


6.8 


10 


1.1 


3.2 


93.6 


0.10 


0.10 


1.5 


12 


3.6 


9 


1.4 


0.9 


7.8 


4.8 


8.8 


0.09 


0.15 


1.7 


run93 


6.5 


2.8 


4 


1.1 


3.2 


47.4 


0.20 


0.20 


1.7 


12 


1.5 


3 


1.4 


1.1 


8.4 


4.8 


4.0 


0.18 


0.30 


1.6 


run94 


6.5 


4.3 


6 


1.1 


3.2 


55.8 


0.20 


0.20 


1.5 


12 


2.2 


5 


1.5 


0.7 


8.3 


4.9 


4.5 


0.18 


0.28 


1.6 


run95 


6.5 


5.7 


8 


1.1 


3.2 


88.0 


0.20 


0.20 


1.7 


12 


3.0 


7 


1.6 


1.3 


8.1 


4.8 


4.7 


0.18 


0.27 


1.6 
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Name Initial Final 





Wo 


M 


N 


Tc 




Pc 


h 


h,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 
hl,ohs 


Pc 


h 


h,c 


c 


run96 


6.5 


7.1 


10 


1.1 


3.2 


101.4 


0.20 


0.20 


1.5 


12 


3.8 


9 


1.6 


1.3 


8.0 


4.8 


6.2 


0.18 


0.26 


1.7 


run97 


7 


2.5 


4 


0.9 


3.2 


61.0 


0.00 


0.00 


1.5 


11 


1.4 


3 


0.3 


0.9 


8.2 


4.6 


249.4 


0.00 


0.00 


2.3 


run98 


7 


3.8 


6 


0.9 


3.3 


73.4 


0.00 


0.00 


1.5 


11 


2.1 


5 


0.3 


0.6 


8.1 


4.8 


320.0 


0.00 


0.00 


2.3 


run99 


7 


5.1 


8 


0.9 


3.3 


110.3 


0.00 


0.00 


1.8 


12 


2.8 


7 


0.4 


0.9 


8.0 


4.5 


281.9 


0.00 


0.00 


2.2 


runlOO 


7 


6.4 


10 


0.9 


3.3 


134.5 


0.00 


0.00 


2.0 


12 


3.5 


9 


0.5 


0.9 


7.9 


4.7 


204.6 


0.00 


0.00 


2.2 


runlOl 


7 


2.6 


4 


0.9 


3.2 


61.1 


0.05 


0.05 


1.5 


12 


1.4 


3 


0.7 


0.3 


8.5 


5.0 


28.3 


0.05 


0.10 


1.9 


runl02 


7 


3.9 


6 


0.9 


3.3 


79.7 


0.05 


0.05 


1.6 


12 


2.1 


5 


0.8 


0.3 


8.4 


5.0 


29.1 


0.05 


0.10 


1.9 


runl03 


7 


5.3 


8 


0.9 


3.3 


118.9 


0.05 


0.05 


1.5 


12 


2.8 


7 


0.9 


0.5 


8.2 


5.0 


33.7 


0.05 


0.10 


1.9 


run 104 


7 


6.6 


10 


0.9 


3.3 


138.3 


0.05 


0.05 


1.6 


12 


3.6 


9 


0.9 


0.5 


8.0 


4.9 


31.2 


0.04 


0.09 


1.9 


runl05 


7 


2.7 


4 


0.9 


3.2 


60.8 


0.10 


0.10 


1.5 


12 


1.4 


3 


1.0 


0.8 


8.7 


5.0 


12.0 


0.09 


0.18 


1.8 


runlOe 


7 


4.0 


6 


0.9 


3.3 


81.2 


0.10 


0.10 


1.6 


12 


2.1 


5 


0.9 


0.4 


8.5 


5.1 


31.0 


0.09 


0.18 


1.9 


runl07 


7 


5.4 


8 


0.9 


3.3 


121.3 


0.10 


0.10 


1.9 


12 


2.9 


7 


1.0 


0.4 


8.3 


5.0 


25.0 


0.09 


0.17 


1.8 


runlOS 


7 


6.8 


10 


0.9 


3.3 


141.0 


0.10 


0.10 


1.6 


12 


3.6 


9 


1.1 


0.5 


8.1 


5.0 


23.5 


0.09 


0.16 


1.8 


runl09 


7 


2.8 


4 


0.9 


3.2 


72.1 


0.20 


0.20 


1.5 


12 


1.5 


3 


1.2 


0.5 


8.8 


5.0 


6.2 


0.18 


0.32 


1.7 


runllO 


7 


4.3 


6 


0.9 


3.3 


83.1 


0.20 


0.20 


1.6 


12 


2.2 


5 


1.1 


0.5 


8.7 


5.2 


13.2 


0.18 


0.30 


1.8 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


fb,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 


Pc 


h 


/6,C 


c 


runlll 


7 


5.7 


8 


0.9 


3.3 


133.2 


0.20 


0.20 


1.6 


12 


3.0 


7 


1.3 


0.6 


8.5 


5.1 


10.4 


0.18 


0.29 


1.7 


runll2 


7 


7.1 


10 


0.9 


3.3 


152.5 


0.20 


0.20 


1.6 


12 


3.8 


9 


1.3 


0.8 


8.3 


5.0 


10.7 


0.18 


0.27 


1.8 


runll3 


7.5 


2.5 


4 


0.7 


3.3 


103.6 


0.00 


0.00 


1.6 


9 


1.4 


3 


0.3 


0.4 


8.3 


4.9 


341.4 


0.00 


0.00 


2.3 


runll4 


7.5 


3.8 


6 


0.7 


3.4 


119.9 


0.00 


0.00 


1.6 


10 


2.1 


5 


0.3 


0.6 


8.2 


4.8 


455.5 


0.00 


0.00 


2.3 


runll5 


7.5 


5.1 


8 


0.7 


3.4 


181.1 


0.00 


0.00 


1.6 


10 


2.8 


7 


0.4 


0.5 


8.1 


4.8 


359.2 


0.00 


0.00 


2.2 


runlie 


7.5 


6.4 


10 


0.7 


3.4 


220.4 


0.00 


0.00 


2.0 


10 


3.5 


9 


0.4 


0.2 


8.0 


4.9 


394.2 


0.00 


0.00 


2.3 


runllT 


7.5 


2.6 


4 


0.7 


3.3 


103.6 


0.05 


0.05 


1.6 


11 


1.4 


3 


0.9 


0.4 


8.7 


4.9 


14.2 


0.05 


0.10 


1.8 


runllS 


7.5 


3.9 


6 


0.7 


3.4 


130.7 


0.05 


0.05 


1.6 


10 


2.1 


5 


0.7 


0.8 


8.6 


4.8 


59.1 


0.05 


0.10 


2.0 


runll9 


7.5 


5.3 


8 


0.7 


3.4 


195.8 


0.05 


0.05 


1.8 


10 


2.8 


7 


0.7 


0.6 


8.4 


5.0 


83.9 


0.05 


0.10 


2.0 


run 120 


7.5 


6.6 


10 


0.7 


3.4 


225.8 


0.05 


0.05 


2.0 


10 


3.6 


9 


0.7 


0.5 


8.2 


5.0 


76.1 


0.05 


0.09 


2.0 


runl21 


7.5 


2.7 


4 


0.7 


3.3 


103.0 


0.10 


0.10 


1.8 


12 


1.4 


3 


0.9 


0.4 


9.0 


5.3 


13.5 


0.09 


0.18 


1.8 


runl22 


7.5 


4.0 


6 


0.7 


3.3 


133.1 


0.10 


0.10 


1.6 


11 


2.1 


5 


0.8 


0.4 


8.8 


5.2 


39.5 


0.09 


0.18 


1.9 


runl23 


7.5 


5.4 


8 


0.7 


3.4 


199.3 


0.10 


0.10 


1.6 


11 


2.9 


7 


0.9 


0.8 


8.5 


5.0 


39.7 


0.09 


0.17 


1.9 


run 124 


7.5 


6.8 


10 


0.7 


3.4 


230.1 


0.10 


0.10 


1.8 


11 


3.6 


9 


0.8 


0.5 


8.4 


5.0 


58.5 


0.09 


0.17 


2.0 


runl25 


7.5 


2.8 


4 


0.7 


3.3 


121.9 


0.20 


0.20 


1.6 


12 


1.4 


3 


1.0 


0.6 


9.2 


5.4 


9.8 


0.18 


0.32 


1.7 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


h,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 
hl,ohs 


Pc 


h 


h,c 


c 


runl26 


7.5 


4.3 


6 


0.7 


3.4 


135.8 


0.20 


0.20 


1.7 


12 


2.2 


5 


0.9 


0.7 


9.1 


5.2 


24.0 


0.18 


0.32 


1.8 


runl27 


7.5 


5.7 


8 


0.7 


3.4 


218.7 


0.20 


0.20 


1.7 


12 


3.0 


7 


1.0 


0.5 


8.9 


5.3 


23.6 


0.18 


0.30 


1.8 


runl28 


7.5 


7.1 


10 


0.7 


3.4 


248.9 


0.20 


0.20 


1.7 


12 


3.8 


9 


1.0 


0.3 


8.8 


5.4 


29.8 


0.18 


0.29 


1.9 


runl29 


4 


3.0 


4 


1.6 


3.3 


14.6 


0.30 


0.30 


1.3 


12 


1.6 


3 


2.4 


2.3 


7.7 


4.4 


0.7 


0.28 


0.36 


1.4 


runlSO 


4 


4.5 


6 


1.6 


3.3 


19.7 


0.30 


0.30 


1.3 


12 


2.4 


5 


2.5 


2.6 


7.5 


4.4 


0.9 


0.27 


0.35 


1.5 


runl31 


4 


6.0 


8 


1.6 


3.3 


28.8 


0.30 


0.30 


1.4 


12 


3.2 


7 


2.6 


2.5 


7.4 


4.5 


1.1 


0.27 


0.34 


1.4 


runl32 


4 


7.5 


10 


1.6 


3.3 


34.7 


0.30 


0.30 


1.7 


12 


4.0 


9 


2.8 


3.1 


7.3 


4.5 


1.1 


0.27 


0.33 


1.5 


runl33 


5 


3.0 


4 


1.5 


3.3 


17.0 


0.30 


0.30 


1.0 


12 


1.6 


3 


2.3 


2.1 


7.8 


4.5 


0.7 


0.27 


0.36 


1.4 


runl34 


5 


4.5 


6 


1.5 


3.3 


22.8 


0.30 


0.30 


1.6 


12 


2.4 


5 


2.4 


2.4 


7.6 


4.5 


1.0 


0.27 


0.36 


1.5 


runl35 


5 


6.0 


8 


1.5 


3.3 


33.7 


0.30 


0.30 


1.2 


12 


3.2 


7 


2.6 


2.4 


7.5 


4.5 


1.1 


0.27 


0.34 


1.5 


runl36 


5 


7.5 


10 


1.5 


3.3 


39.8 


0.30 


0.30 


1.8 


12 


4.0 


9 


2.6 


2.7 


7.3 


4.4 


1.5 


0.27 


0.34 


1.5 


runl37 


6 


3.0 


4 


1.4 


3.2 


20.6 


0.30 


0.30 


1.1 


12 


1.6 


3 


2.2 


2.1 


7.9 


4.4 


1.0 


0.28 


0.37 


1.4 


runl38 


6 


4.5 


6 


1.4 


3.3 


27.1 


0.30 


0.30 


1.4 


12 


2.4 


5 


2.3 


1.9 


7.7 


4.5 


1.1 


0.27 


0.36 


1.5 


runl39 


6 


6.0 


8 


1.4 


3.3 


40.4 


0.30 


0.30 


1.5 


12 


3.2 


7 


2.5 


2.2 


7.6 


4.5 


1.3 


0.27 


0.34 


1.5 


runl40 


6 


7.5 


10 


1.4 


3.3 


47.6 


0.30 


0.30 


1.7 


12 


4.0 


9 


2.5 


2.5 


7.5 


4.5 


1.6 


0.27 


0.34 


1.5 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


fb,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 


Pc 


h 


/6,C 


c 


runl41 


7 


3.0 


4 


1.3 


3.2 


25.8 


0.30 


0.30 


1.2 


12 


1.5 


3 


2.0 


1.5 


8.1 


4.6 


1.1 


0.28 


0.37 


1.5 


runl42 


7 


4.5 


6 


1.3 


3.2 


33.4 


0.30 


0.30 


1.2 


12 


2.4 


5 


2.2 


1.7 


7.9 


4.5 


1.5 


0.27 


0.36 


1.5 


runl43 


7 


6.0 


8 


1.3 


3.3 


50.3 


0.30 


0.30 


1.8 


12 


3.2 


7 


2.1 


1.7 


7.7 


4.6 


2.1 


0.27 


0.36 


1.5 


runl44 


7 


7.5 


10 


1.3 


3.2 


58.8 


0.30 


0.30 


1.8 


12 


4.0 


9 


2.4 


2.0 


7.6 


4.6 


1.8 


0.27 


0.34 


1.5 


runl45 


8 


3.0 


4 


1.2 


3.2 


34.0 


0.30 


0.30 


1.3 


12 


1.5 


3 


1.8 


1.3 


8.3 


4.7 


1.6 


0.28 


0.39 


1.5 


runl46 


8 


4.5 


6 


1.2 


3.2 


43.2 


0.30 


0.30 


1.3 


12 


2.4 


5 


1.9 


1.7 


8.1 


4.7 


2.1 


0.27 


0.37 


1.5 


runl47 


8 


6.0 


8 


1.2 


3.2 


66.3 


0.30 


0.30 


1.5 


12 


3.2 


7 


2.0 


1.7 


7.9 


4.7 


2.5 


0.27 


0.36 


1.6 


runl48 


8 


7.5 


10 


1.2 


3.2 


76.0 


0.30 


0.30 


1.7 


12 


4.0 


9 


2.1 


1.6 


7.8 


4.7 


2.7 


0.27 


0.35 


1.6 


runl49 


6.5 


3.0 


4 


1.1 


3.2 


47.8 


0.30 


0.30 


1.5 


12 


1.5 


3 


1.5 


1.1 


8.6 


4.8 


2.8 


0.28 


0.40 


1.6 


run 150 


6.5 


4.5 


6 


1.1 


3.2 


59.5 


0.30 


0.30 


1.7 


12 


2.3 


5 


1.8 


1.3 


8.4 


4.8 


2.7 


0.27 


0.37 


1.6 


runl51 


6.5 


6.0 


8 


1.1 


3.2 


92.1 


0.30 


0.30 


1.5 


12 


3.2 


7 


1.7 


1.3 


8.2 


4.9 


4.3 


0.27 


0.37 


1.6 


runl52 


6.5 


7.5 


10 


1.1 


3.2 


104.6 


0.30 


0.30 


1.7 


12 


4.0 


9 


1.7 


1.4 


8.1 


4.8 


5.1 


0.27 


0.37 


1.6 


runl53 


7 


3.0 


4 


0.9 


3.2 


72.5 


0.30 


0.30 


1.5 


12 


1.5 


3 


1.4 


0.8 


9.0 


5.1 


3.9 


0.28 


0.42 


1.6 


run 154 


7 


4.5 


6 


0.9 


3.3 


88.3 


0.30 


0.30 


1.8 


12 


2.3 


5 


1.4 


0.9 


8.8 


5.1 


6.2 


0.27 


0.40 


1.7 


runl55 


7 


6.0 


8 


0.9 


3.3 


139.3 


0.30 


0.30 


1.5 


12 


3.1 


7 


1.4 


0.9 


8.6 


5.0 


7.2 


0.27 


0.38 


1.7 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


fb,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 


Pc 


h 


/6,C 


c 


runl56 


7 


7.5 


10 


0.9 


3.3 


157.0 


0.30 


0.30 


1.7 


12 


3.9 


9 


1.6 


1.0 


8.5 


5.0 


6.7 


0.27 


0.37 


1.7 


runl57 


7.5 


3.0 


4 


0.7 


3.3 


122.3 


0.30 


0.30 


1.6 


12 


1.5 


3 


1.1 


0.4 


9.4 


5.4 


10.0 


0.28 


0.43 


1.7 


runl58 


7.5 


4.5 


6 


0.7 


3.3 


143.3 


0.30 


0.30 


1.6 


12 


2.3 


5 


1.1 


0.5 


9.3 


5.4 


11.6 


0.27 


0.42 


1.8 


runl59 


7.5 


6.0 


8 


0.7 


3.4 


229.1 


0.30 


0.30 


1.6 


12 


3.1 


7 


1.2 


0.6 


9.1 


5.4 


13.3 


0.27 


0.40 


1.8 


runieO 


7.5 


7.5 


10 


0.7 


3.4 


255.8 


0.30 


0.30 


1.7 


12 


3.9 


9 


1.2 


0.7 


9.0 


5.4 


15.5 


0.27 


0.39 


1.8 


runiei 


4 


2.6 


4 


1.2 


2.5 


29.1 


0.05 


0.05 


1.3 


12 


1.4 


3 


0.9 


0.6 


6.1 


3.3 


12.1 


0.05 


0.10 


1.8 


runl62 


4 


3.9 


6 


1.2 


2.5 


41.0 


0.05 


0.05 


1.4 


12 


2.1 


5 


1.2 


0.8 


5.8 


3.4 


8.5 


0.05 


0.08 


1.7 


runl63 


4 


5.3 


8 


1.2 


2.5 


58.5 


0.05 


0.05 


1.7 


12 


2.8 


7 


1.4 


1.1 


5.7 


3.3 


7.2 


0.05 


0.08 


1.7 


runl64 


4 


6.6 


10 


1.2 


2.5 


71.7 


0.05 


0.05 


1.8 


12 


3.6 


9 


1.5 


1.2 


5.6 


3.3 


6.9 


0.04 


0.08 


1.7 


runl65 


4 


2.7 


4 


1.2 


2.5 


29.3 


0.10 


0.10 


1.4 


12 


1.4 


3 


1.0 


0.6 


6.2 


3.3 


11.2 


0.09 


0.17 


1.8 


runl66 


4 


4.0 


6 


1.2 


2.5 


41.9 


0.10 


0.10 


1.7 


12 


2.2 


5 


1.2 


1.0 


6.0 


3.4 


7.5 


0.09 


0.16 


1.7 


runl67 


4 


5.4 


8 


1.2 


2.5 


60.0 


0.10 


0.10 


2.0 


12 


2.9 


7 


1.5 


1.0 


5.8 


3.3 


6.0 


0.09 


0.15 


1.7 


runl68 


4 


6.8 


10 


1.2 


2.5 


73.2 


0.10 


0.10 


1.5 


12 


3.6 


9 


1.5 


0.9 


5.7 


3.3 


6.2 


0.09 


0.14 


1.7 


runl69 


4 


2.8 


4 


1.2 


2.5 


33.9 


0.20 


0.20 


1.0 


12 


1.5 


3 


1.2 


0.8 


6.4 


3.4 


4.9 


0.18 


0.29 


1.7 


runl70 


4 


4.3 


6 


1.2 


2.5 


43.3 


0.20 


0.20 


1.6 


12 


2.3 


5 


1.4 


1.1 


6.1 


3.4 


5.6 


0.18 


0.29 


1.7 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


h,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 


Pc 


h 


/6,C 


c 


runlTl 


4 


5.7 


8 


1.2 


2.5 


65.4 


0.20 


0.20 


1.7 


12 


3.0 


7 


1.6 


1.4 


6.0 


3.4 


4.8 


0.18 


0.27 


1.7 


runl72 


4 


7.1 


10 


1.2 


2.5 


79.2 


0.20 


0.20 


1.8 


12 


3.8 


9 


1.6 


1.3 


5.8 


3.4 


5.3 


0.18 


0.26 


1.7 


runl73 


4 


3.0 


4 


1.2 


2.5 


34.6 


0.30 


0.30 


0.9 


12 


1.5 


3 


1.4 


0.9 


6.6 


3.5 


3.8 


0.27 


0.40 


1.6 


runl74 


4 


4.5 


6 


1.2 


2.5 


46.6 


0.30 


0.30 


1.3 


12 


2.4 


5 


1.5 


1.3 


6.3 


3.4 


4.3 


0.27 


0.38 


1.6 


runl75 


4 


6.0 


8 


1.2 


2.5 


68.3 


0.30 


0.30 


1.3 


12 


3.2 


7 


1.6 


1.3 


6.1 


3.3 


4.7 


0.27 


0.37 


1.7 


runl76 


4 


7.5 


10 


1.2 


2.5 


82.2 


0.30 


0.30 


1.3 


12 


4.0 


9 


1.8 


1.5 


6.0 


3.4 


4.5 


0.27 


0.35 


1.7 


runl77 


5 


2.6 


4 


1.1 


2.4 


41.1 


0.05 


0.05 


1.3 


12 


1.4 


3 


0.6 


0.6 


6.3 


3.4 


36.4 


0.05 


0.10 


1.9 


runl78 


5 


3.9 


6 


1.1 


2.5 


56.8 


0.05 


0.05 


1.2 


12 


2.1 


5 


1.0 


0.6 


6.1 


3.5 


14.9 


0.05 


0.09 


1.8 


runl79 


5 


5.3 


8 


1.1 


2.5 


82.0 


0.05 


0.05 


1.5 


12 


2.8 


7 


1.2 


0.8 


5.9 


3.4 


10.9 


0.04 


0.08 


1.8 


run 180 


5 


6.6 


10 


1.1 


2.5 


98.6 


0.05 


0.05 


1.8 


12 


3.6 


9 


1.2 


1.0 


5.8 


3.4 


12.2 


0.04 


0.08 


1.8 


runlBl 


5 


2.7 


4 


1.1 


2.4 


41.2 


0.10 


0.10 


1.3 


12 


1.4 


3 


0.8 


0.4 


6.5 


3.6 


20.1 


0.09 


0.18 


1.9 


runl82 


5 


4.0 


6 


1.1 


2.5 


58.0 


0.10 


0.10 


1.4 


12 


2.2 


5 


1.0 


0.6 


6.2 


3.4 


12.8 


0.09 


0.17 


1.8 


runlBS 


5 


5.4 


8 


1.1 


2.5 


84.0 


0.10 


0.10 


1.2 


12 


2.9 


7 


1.2 


1.0 


6.0 


3.5 


9.9 


0.09 


0.15 


1.8 


run 184 


5 


6.8 


10 


1.1 


2.5 


100.7 


0.10 


0.10 


1.8 


12 


3.6 


9 


1.3 


1.2 


5.9 


3.4 


10.3 


0.09 


0.15 


1.8 


runl85 


5 


2.8 


4 


1.1 


2.4 


48.2 


0.20 


0.20 


1.2 


12 


1.5 


3 


1.0 


0.5 


6.7 


3.6 


9.9 


0.18 


0.31 


1.8 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 


Tc 




Pc 


h 


h,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 
hl,ohs 


Pc 


h 


h,c 


c 


runl86 


5 


4.3 


6 


1.1 


2.5 


59.8 


0.20 


0.20 


1.7 


12 


2.3 


5 


1.2 


0.9 


6.4 


3.5 


8.8 


0.18 


0.29 


1.7 


runl87 


5 


5.7 


8 


1.1 


2.5 


91.8 


0.20 


0.20 


1.4 


12 


3.0 


7 


1.4 


0.9 


6.2 


3.5 


6.7 


0.18 


0.27 


1.7 


runl88 


5 


7.1 


10 


1.1 


2.5 


109.0 


0.20 


0.20 


1.9 


12 


3.8 


9 


1.4 


0.9 


6.1 


3.4 


8.6 


0.18 


0.27 


1.8 


runl89 


5 


3.0 


4 


1.1 


2.4 


48.9 


0.30 


0.30 


1.1 


12 


1.5 


3 


1.2 


0.9 


6.8 


3.5 


5.9 


0.27 


0.41 


1.7 


runl90 


5 


4.5 


6 


1.1 


2.4 


64.2 


0.30 


0.30 


1.6 


12 


2.3 


5 


1.4 


1.2 


6.6 


3.5 


5.9 


0.27 


0.39 


1.7 


runl91 


5 


6.0 


8 


1.1 


2.5 


95.9 


0.30 


0.30 


1.2 


12 


3.2 


7 


1.5 


1.1 


6.3 


3.4 


5.8 


0.27 


0.37 


1.7 


runl92 


5 


7.5 


10 


1.1 


2.5 


112.8 


0.30 


0.30 


1.9 


12 


4.0 


9 


1.5 


0.9 


6.2 


3.4 


6.8 


0.27 


0.36 


1.7 


runl93 


6 


2.6 


4 


0.9 


2.4 


67.8 


0.05 


0.05 


1.4 


11 


1.4 


3 


0.5 


0.4 


6.7 


3.5 


79.2 


0.05 


0.11 


2.1 


run 194 


6 


3.9 


6 


0.9 


2.4 


91.3 


0.05 


0.05 


1.3 


12 


2.1 


5 


0.7 


0.6 


6.5 


3.5 


53.2 


0.05 


0.10 


2.0 


runl95 


6 


5.3 


8 


0.9 


2.4 


134.1 


0.05 


0.05 


1.5 


12 


2.8 


7 


0.7 


0.4 


6.3 


3.6 


81.8 


0.04 


0.10 


2.0 


runl96 


6 


6.6 


10 


0.9 


2.4 


157.7 


0.05 


0.05 


1.7 


12 


3.6 


9 


0.8 


0.5 


6.2 


3.6 


56.5 


0.04 


0.09 


2.0 


runl97 


6 


2.7 


4 


0.9 


2.4 


67.8 


0.10 


0.10 


1.3 


12 


1.4 


3 


0.6 


0.4 


6.9 


3.9 


49.9 


0.09 


0.18 


2.0 


runl98 


6 


4.0 


6 


0.9 


2.4 


93.2 


0.10 


0.10 


1.3 


12 


2.2 


5 


0.8 


0.4 


6.6 


3.7 


36.1 


0.09 


0.18 


1.9 


runl99 


6 


5.4 


8 


0.9 


2.4 


137.2 


0.10 


0.10 


1.7 


12 


2.9 


7 


0.8 


0.5 


6.4 


3.7 


35.0 


0.09 


0.17 


1.9 


run200 


6 


6.8 


10 


0.9 


2.4 


161.0 


0.10 


0.10 


1.7 


12 


3.6 


9 


1.0 


0.6 


6.3 


3.6 


27.8 


0.09 


0.16 


1.9 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


h,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 
hl,ohs 


Pc 


h 


h,c 


c 


run201 


6 


2.8 


4 


0.9 


2.4 


79.8 


0.20 


0.20 


1.3 


12 


1.5 


3 


0.7 


0.4 


7.1 


3.9 


33.7 


0.18 


0.33 


1.9 


run202 


6 


4.3 


6 


0.9 


2.4 


95.7 


0.20 


0.20 


1.2 


12 


2.2 


5 


0.9 


0.6 


6.8 


3.8 


18.3 


0.18 


0.31 


1.8 


run203 


6 


5.7 


8 


0.9 


2.4 


150.3 


0.20 


0.20 


1.6 


12 


3.0 


7 


1.0 


0.6 


6.6 


3.7 


21.7 


0.18 


0.29 


1.8 


run204 


6 


7.1 


10 


0.9 


2.4 


174.4 


0.20 


0.20 


1.7 


12 


3.8 


9 


1.2 


0.8 


6.4 


3.6 


16.5 


0.18 


0.28 


1.8 


run205 


6 


3.0 


4 


0.9 


2.4 


80.5 


0.30 


0.30 


1.3 


12 


1.5 


3 


1.0 


0.7 


7.3 


3.8 


10.2 


0.27 


0.42 


1.8 


run206 


6 


4.5 


6 


0.9 


2.4 


102.4 


0.30 


0.30 


1.6 


12 


2.3 


5 


1.1 


0.4 


7.0 


3.7 


11.9 


0.27 


0.40 


1.8 


run207 


6 


6.0 


8 


0.9 


2.4 


157.1 


0.30 


0.30 


1.3 


12 


3.1 


7 


1.2 


1.1 


6.7 


3.6 


10.8 


0.27 


0.39 


1.8 


run208 


6 


7.5 


10 


0.9 


2.4 


180.2 


0.30 


0.30 


1.8 


12 


3.9 


9 


1.3 


0.9 


6.6 


3.7 


11.9 


0.27 


0.37 


1.8 


run209 


7 


2.6 


4 


0.7 


2.4 


144.9 


0.05 


0.05 


1.5 


9 


1.4 


3 


0.7 


0.4 


7.0 


3.8 


34.5 


0.05 


0.11 


1.9 


run210 


7 


3.9 


6 


0.7 


2.4 


189.0 


0.05 


0.05 


1.5 


10 


2.1 


5 


0.6 


0.5 


6.8 


3.7 


69.0 


0.05 


0.10 


2.0 


run211 


7 


5.3 


8 


0.7 


2.5 


281.9 


0.05 


0.05 


1.5 


9 


2.9 


7 


0.6 


0.7 


6.5 


3.5 


91.5 


0.05 


0.10 


2.0 


run212 


7 


6.6 


10 


0.7 


2.5 


327.7 


0.05 


0.05 


1.7 


9 


3.6 


9 


0.7 


0.6 


6.4 


3.6 


88.1 


0.04 


0.09 


2.0 


run213 


7 


2.7 


4 


0.7 


2.4 


144.2 


0.10 


0.10 


1.5 


10 


1.4 


3 


0.8 


0.4 


7.2 


4.0 


21.7 


0.09 


0.19 


1.9 


run214 


7 


4.0 


6 


0.7 


2.4 


192.6 


0.10 


0.10 


1.5 


9 


2.2 


5 


0.8 


0.3 


6.9 


3.8 


38.9 


0.09 


0.18 


1.9 


run215 


7 


5.4 


8 


0.7 


2.5 


287.6 


0.10 


0.10 


1.6 


10 


2.9 


7 


0.7 


0.3 


6.7 


3.8 


76.3 


0.09 


0.18 


2.0 
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Table 1 - Continued 



Name Initial Final 





Wo 


M 


N 






Pc 


h 


fb,c 


c 


t 


M 


N 




Ar,cut 
c,obs 




„cut 


Pc 


h 


/6,C 


c 


run216 


7 


6.8 


10 


0.7 


2.5 


334.2 


0.10 


0.10 


1.5 


10 


3.7 


9 


0.7 


0.7 


6.6 


3.7 


69.2 


0.09 


0.17 


2.0 


run217 


7 


2.8 


4 


0.7 


2.4 


171.0 


0.20 


0.20 


1.5 


11 


1.5 


3 


0.9 


0.7 


7.5 


3.9 


18.6 


0.18 


0.32 


1.8 


run218 


7 


4.3 


6 


0.7 


2.4 


197.0 


0.20 


0.20 


1.5 


10 


2.3 


5 


0.9 


0.4 


7.2 


4.0 


26.9 


0.18 


0.31 


1.9 


run219 


7 


5.7 


8 


0.7 


2.5 


315.7 


0.20 


0.20 


1.6 


11 


3.0 


7 


0.8 


0.6 


7.0 


3.8 


42.3 


0.18 


0.31 


1.9 


run220 


7 


7.1 


10 


0.7 


2.5 


361.5 


0.20 


0.20 


1.8 


11 


3.8 


9 


0.8 


0.5 


6.9 


3.9 


51.0 


0.18 


0.29 


2.0 


run221 


7 


3.0 


4 


0.7 


2.4 


171.8 


0.30 


0.30 


1.5 


12 


1.5 


3 


1.0 


0.6 


7.9 


4.3 


12.4 


0.27 


0.42 


1.8 


run222 


7 


4.5 


6 


0.7 


2.4 


209.2 


0.30 


0.30 


1.6 


12 


2.3 


5 


0.9 


0.5 


7.6 


4.1 


23.9 


0.27 


0.42 


1.9 


run223 


7 


6.0 


8 


0.7 


2.5 


330.2 


0.30 


0.30 


1.5 


11 


3.1 


7 


0.9 


0.6 


7.3 


4.0 


26.5 


0.27 


0.41 


1.9 


run224 


7 


7.5 


10 


0.7 


2.4 


372.1 


0.30 


0.30 


1.9 


11 


3.9 


9 


1.0 


0.6 


7.1 


3.9 


31.2 


0.27 


0.39 


1.9 
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